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Stellar evolution tracks and isochrones are key inputs for a wide range of 
astrophysical studies; in particular, they are essential to the interpretation of 
photometric and spectroscopic observations of resolved and unresolved stellar 
populations. We have made available to the astrophysical community a large, 
homogenous database of up-to-date stellar tracks and isochrones, and a set of 
programs useful in population synthesis studies. In this paper we first summa- 



rize the main properties o f our stella r model database 



duced in iPietrinferni et al.l (120041 ) and iPietrinferni et al.l (120061 ) . We then discuss 



'BaSTI) already intro- 



an important update of the database, i.e., the extension of all stellar models 
and isochrones until the end of the thermal pulses along the Asymptotic Giant 
Branch. This extension of the library is particularly relevant for stellar popula- 
tion analyses in the near-infrared, or longer wavelengths, where the contribution 
to the integrated photometric properties by cool and bright Asymptotic Giant 
Branch stars is significant. A few comparisons with empirical data are also pre- 
sentend and briefly discussed. We then present three web-tools that allow an 
interactive access to the database, and make possible to compute user-specified 
evolutionary tracks, isochrones, stellar luminosity functions, plus synthetic Color- 
Magnitude-Diagrams and integrated magnitudes for arbitrary Star Formation 
Histories. All these web tools are available at the BaSTI database official site: 
http : / /www . oa-teramo . inaf . it/BASTlj 



Subject headings: astronomical data bases: miscellaneous - galaxies: stellar con- 
tent - stars: general - stars: AGB and post-AGB 



1. Introduction 



The interpretation of photometric and spectroscopic observations of resolved and unre- 
solved stellar populations is nowadays a fundamental tool to investigate the formation and 
evolution of galaxies. Large grids of stellar evolution models and isochrones are a neces- 
sary ingredient in this kind of analyses, together with appropriate tools to predict synthetic 
Color-Magnitude-Diagrams (CMDs) - hence star counts along the various observed CMD 
branches - integrated magnitudes and spec tra of stellar popu l ations with an arbitrary Star 



Forma tion History (SFH). To this purpose IPietrinferni et al.l (120041 ) and IPietrinferni et al. 



( 120061 ) have published a huge and homogeneous database of stellar evolution models, that 
covers the relevant chemical composition range of stellar populations in galaxies of various 
morphological types, and allows one to choose among different treatments of core convection 
and stellar mass loss employed in the model calculations. These models are available on the 
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web as downloadable files, from the BaSTI database - an acronym for "a Bag of Stellar 
Tracks and Isochrones" - at |http: //www, oa-teramo . inaf . it/BASTl] , together with three 
web-tools that allow user-friendly manipulations of the stellar evolution library. Thanks to 
these tools one can compute user-specified isochrones, interpolate among evolutionary tracks 
to obtain the evolution for a mass value not contained in the grid, determine the differential 
and cumulative luminosity functions from a set of isochrones and, finally, compute synthetic 
stellar populations with an arbitrary SFH and determine their integrated magnitudes and 
colors. 

This paper describes comprehensively the BaSTI database, and in particular how to 
run the various web-tools, the necessary inputs and their outputs. Next section describes 
the stella r evolution library , with parti cular emphasis on the m ost recent extensions not in- 
cluded in iPietrinferni et al.l (120041 ) and iPietrinferni et al.l (120061 ); the following three sections 
describe the individual web-tools that allow compute user specified isochrones and models, 
luminosity functions and synthetic CMDs. Final conclusions and an appendix with some 
more technical details about the implementation of the web-tools close the paper. 



2. Stellar evolution models and isochrones 



We have computed stellar evolution models for 11 different metallicities, namely Z = 
0.0001, 0.0003, 0.0006, 0.001, 0.002, 0.004, 0.008, 0.01, 0.0198 0.03 and 0.04, assumin g 
two different heavy element distributions: the scaled-s olar one by iGreyesse fc Noelsl (119931 ). 
and the a-enhanced one with < [a/Fe] >= 0.40 by ISalaris fc Weissl (119981 ). As for the 
initi al He-mass fr a ction s, we employed a value Y = 0.245 for the cosmological He (see, 



e.g 



Cassisi et al.l (120031 ) ) . To reproduce the solar initial He-abundance obtained from the 
calibration of the solar model, we assume an Helium enrichment law equal to AY/AZ ~ 1.4. 
For each initial chemical composition we have computed stellar evolutionary tracks with mass 
in the range 0.5 < M/ M^ < 10 and a fi ne mass spacing. Mass loss from the stellar surface is 
accounted for using the iReimerd (119751 ) law and two values of the free parameter rj, namely 
77=0.2 and 0.4. For stars that develop convective cores during the central H-burning phase 
we have computed models with and without overshoot from the Schwarzschild boundary of 
the central convective regions. 

All evolutions, with the exception of the least massive stars whose central H-burning 
timescale is longer than the Hubble time, have been computed from the Pre-Main Sequence 
phase until to the C-ignition, or until the first thermal pulse along the Asymptotic Giant 
Branch (AG B), depending on t he in itial s tellar mass. Full details of these calculations can 
be found in IPietrinferni et al.l (120041 ) and IPietrinferni et al.l ( 120061 ). Very recently we have 
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extended the evolution of low- and intermediate-mass objects through the entire Thermal 
Pulse (TP) phase, until the end of the AGB stage. To this purpose we have used the 
synthetic AGB technique ( Iben Sz Truran ( 19781 )). Our simplified treatment of the TP- AGB 
phase reproduces satisfactorily several integrated properties (see below) of stellar populations 
in the near-IR bands, which are greatly affected by the presence of TP stars. Given that 
the purpose of this library of stellar models is to provide a reliable tool to be employed in 
stellar population synthesis studies, we feel confident that the following simplified treatment 
of the AGB evolution is adequate for our purposes. For each stellar model of a given initial 
chemical composition and mass, we started the synthetic AGB evolution at the beginning 
of the TP phase, where the stellar evolution calculations were terminated. The first model 
of the synthetic TP-AGB evolution is characterized by the total mass M, Carbon-Oxygen 
core mass Mco, luminosity L, effective temperature T e ff and surface chemical composition 
(Xf, Yj, Zj) of the last fully evolutionary model. 

The TP-AGB phas e is then followed by increasing (after a given timestep dt) Mco and L 
according to Eqs. 5-7 in lWagenhuber fc Groenewegenl (119981 ) which contain a term mimicking 
the effect of the Hot Bottom Bur ning;, when appropriate. The hydro gen mass fraction in 
the envelope (an input of Eq. 6 of IWagenhuber fc Groenewegenl (119981 )) is approximated as 
1 — (Yf + 0.01) — Zf all along the TP evolution. During a timestep the mass of the envelope 
(M — Mco) is reduced not only by the growt h of Mco, but also b y ma ss loss processes, 
accounted for using the mass loss formulae by IVassiliadis fc Woodl (119931 ) . For any given 
value of M and M co, the effective temperatures are computed using the rela ti onshi ps in 
Wagenhuberl (119961 ). that are plotted in Fig. 8 of I Wagenhuber &: Groenewegenl (119981 ). To 
ensure continuity, the zero points of the equations describing the evolution of L, Mco and 
T e ff have been adjusted to reproduce the corresponding values of the last fully evolutionary 
model, at the beginning of the TP phase. The synthetic evolution is stopped when the 
models have started to evolve off the AGB, at constant luminosity, towards their White 
Dwarf cooling sequence. The full models have been used to compute isochrones for ages 
ranging between 30 Myr and 19 Gyr. We provide separately models and isochrones with and 
without the inclusion of the TP-AGB phase, in case users wish to use their own description of 
the TP-AGB evolution, either employing fully evolutionary calculations, or synthetic AGB 
models with different choices about, i.e., mass loss and/or core mass-luminosity relationship. 
To this purpose we provide - as mentioned below - tables with He and CO core mass as well 
as the luminosity and effective temperature at the beginning of the TP-AGB phase. 



The reader is referred to iPietrinferni et al.l (120041 ) for a detailed discussion about the 
normalization of the individual evolutionary tracks (the concept of key-point and normalized 
points between two successive key-points) and iso chrone computation. He re we simply men- 
tion two changes with respect to the discussion in IPietrinferni et al.l (120041 ) . The first change 
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is the inclusion of two additional key-points along the Red Giant Branch (RGB) one to mark 
the maximum luminosity during the RGB bump, and another one to mark the minimum 
RGB bump luminosity (key points 6 and 7). A total of 370 normalized points is distributed 
between key point 6 and its predecessor, 30 points are distributed between key-points 6 and 
7, and 400 points are distributed between key-point 7 and the following one, the tip of the 
RGB. Due to the TP-AGB extension, we have added an additional key-point (number 17) 
after the beginning of the TP phase, to mark the termination of the AGB phase, when the 
model evolution turns to the blue. A total of 250 normalized points cover the TP-AGB 
phase, between the key-points 16 and 17. The resulting isochrones have a total of 2250 
points. 

The main characteristics of the stellar model database are summarized in Tables 1 and 
2. Table 1 lists the grid of initial chemical compositions; Table 2 summarizes the available 
sets of models and isochrones. For a given heavy elements mixture (either scaled-solar or a- 
enhanced) we have computed models for 11 Z values (Table 1) using both rj=0.2 and i]=0A 
for the evolution until the TP-AGB phase. For each 77 we have computed models with and 
without overshoot from the central convective core (in the appropriate mass range). The 
extension of the overshooting region \ ov is 0.2H P (H p is the local pressure scale height at 
the Schwarzschi ld convective boundary) decreasing to zero when the convective cores vanish, 
as described in IPietrinferni et al.l (120041 ) . For each choice of the metal mixture, Z, rj and 
convective core extension we specify the number of tracks available, their mass range, the 
number of isochrones available and their age range. Notice that for masses above ~ 2.4M 
the evolution with r/=0.4 is identical to the case of r]=0.2 (because of the negligible amount 
of mass lost during their evolution, according to the Reimers' law). 

Broadband magnitudes and colors of the stellar evolution tracks and isochrones are 
predicted using color- T e ff transformatio ns and bolometric correc tions based on an update d 
set of model atmospheres described in IPietrinferni et al.l (120041 ) and ICassisi et al.l (120041 ) . 
The evolutionary results are available in the photometric filters listed in Table 2. Additional 
filters will be added with time. It is worth to point out that, for the first time, we have 
homogenoeus transformations for both scaled-solar and a— enhanced mixture even for super- 
solar metallicities. At present, the transformations for the ACS HST filters are available only 
for a scaled-solar chemical composition. 

For the TP-AGB section of the models and isochrones broadband colors and magnitudes 
have been compute d by supplementing t he transformations used up to the beginning of the 
TP phase with the IWestera et aD (booj l ones for RGB and AGB stars with T P . ff < 3750 K. 
To ensure continuity, the bolometric corrections and colors of l Westera et al.l (120021 ) have been 
shifted to match the other sets of transformations when T e //=3750 K. Along the TP-AGB 
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phase, when the (J — K) colors r each {J — K)=1.2 m ap;, we used the (J — K) — T e ff and 
(H — K) — T e ff relationships by iBergeat et al.l (120011 ) that are appropriate in the carbon 
star regime. Models and isochrones extended along the TP-AGB phase are at the moment 
transformed only to the UBVRIJHKL system. 

All models described above have been computed in the assumption that a star of a 
given mass and initial chemical composition loses during its evolution a fixed amount of 
mass, determined by the adopted mass loss law. An isochrone for an old population will 
therefore display on the CMD a clumpy Horizontal Branch (HB) phase, made of stars all 
essentially with the same total mass, because of both the constant mass of their progenitor, 
and the constant value of the mass lost along the previous RGB phase. If one wants to 
simulate the extended HBs observed in Galactic globular clusters, a spread in the amount of 
mass lost by the RGB progenitor has to be accounted for. To this purpose, for each chemical 
composition we have also computed additional large sets of core He-burning models, with 
He core mass and envelope chemical profile fixed by a RGB progenitor having an age of 
~ 13 Gyr at the RGB tip, and a range of values of the total stellar mass, to simulate the 
effect of a spread in the mass lost along the RGB. These Horizontal Branch (HB) models 
(~ 30 for each chemical composition) constitute a valuable tool to perform synthetic HB 
modeling, and to investigate pulsational and evolutionary properties of different kinds of 
variable stars. 

Users can directly download from the BaSTI database evolutionary tracks and HB 
models as single files, or as tar gzipped archive files, selecting among the following options 
(11 initial chemical compositions for each choice): scaled-solar, a-enhanced, canonical (no- 
overshoot), non-canonical (overshoot included), r/=0.2, r]=0A, TP-AGB excluded, TP-AGB 
included. We have included in the web page appropriate readme files and provide also tables 
(both in ascii and html format) summarizing the main properties of the theoretical models 
(not including the TP-AGB phase) and providing also the CMD location of the Zero Age 
HB (ZAHB) and the central He exhaustion stages. These two latter sets of tables can be 
found at the pages containing the 77=0.4 models. A large number of precomputed isochrones 
for each of the individual options described before - spanning the full age range allowed by 
the computed models - can be downloaded as tar gzipped archive files. 



The m odels and isochrone s conta ined in BaSTI have been extensively tested in lPietrinferni et al 



(12004) and iPietrinferni et al.l (120061 ). and already employed in a number of investigations 
(e.g., iGallart et al.l (120051 )). We conclude this section showing a few tests for the TP-AGB 
extension we introduced in this paper. We have first tested our isochrones against near-IR 
Surface-Brightness-Fluctuation (SBF) magnitudes, that are very se n sitive to the brightnes s 
and evolutionary timescales of AGB stars (see, e.g.. iLiu et al.l (120001 ). ICantiello et al.l (120031 ). 
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Mouhcine et al.l (12005( 1 . iRaimondo et all (120051 ) ). 



We considered the J and .fT-band SBF magnitudes determined by lGonzalez et al.l (120041 ) 
for a sample of Magellanic Clouds' superclusters. A s upercluster i s obta ined by coadding 



individual clusters belonging to the same SWB class f seelSearle et al. 



of the SWB classes, that are essentially age ranges). IGonzalez et al. 



(198PJ) for the definition 
(12004 ) provide also the 



mean metallicity and age of the objects grouped in each supercluster. The original SBF 
magnitu des were on the 2MASS syst em and ha v e bee n transformed to the Johnson system 
following liessell k Brettl <jl988h and ICarpenterl (boOlh . 



Figures [T]and [2] display the observed SBF absolute magnitudes in K and J (IGonzalez et al. 

d2004h assume (m - M) = 18.50 ± 0.13 for the Large Magellanic Cloud (LMC) and 
(m — M)q = 18.99 ± 0.05 for the SMC) of 6 superclusters - whose age range is covered 
by our models - as a function of age. Different symbols correspond to different mean metal- 
licities; the three youngest superclusters share all the same metallicity. Horiz ontal error bars 



corres pond to the age spread around the mean supercluster ages, as given by lGonzalez et al. 
(h004h . 



We have compared these data with the SBF magnitudes obtained from our isochrones, 
for the appropriate mean metallicities of the individual superclusters. Notice that for each 
individual Z considered, we have plotted the theoretical SBF magnitudes only in the age 
range spanned by the one (or more) supercluster with that mean value of Z. We employed 
the BaSTI models with rj—0.2, including overshooting in the appropriate age range. For 
the lowest metallicity (oldest) supercluster we display also the results obtained with rj=0A. 
Theory is able to fit the data points within the observational errors, the only exception being 
the J-band SBF of the oldest supercluster. 

Figure B] compares t he dereddened integrated (J — K) and (H — K) colors provided 
by IGonzalez et al.l (120041 ) for the 6 superclusters - transformed to the Johnson system as 
done for the SBF - with theoretical predictions from the same models employed for the 
SBF analysis. Theoretical colors appear consistent with the observational counterpart; the 
mean difference between theory and observations is equal to ~0.03 mag in both (J — K) 
and (H — K). This value can be considered a reasonable estimate of the accuracy of our 
predicted integrated near-IR colors. Finally, Fig. H] compares the combined K — (J — K) 
CMDs of AGB stars in a sample of LMC clusters, with theoretical isochrones. The individual 
stars have been coadded acc ording to the SWB class of the parent clusters. Reddenings are 
from lFrogel fc Cohenl (119821 ) and we used an LMC distance modulus (m — M) =18.50. The 
metallicities and ages of t he theoretical is o chron es have been selected on the basis of the 
SWB class following again IGonzalez et al.l (120041 ) . One can notice a good overlap between 
the theoretical AGB sequences and the observed near-IR CMDs. 
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3. Isochrone Maker and Track interpolation program 

We have realized a dedicated WEB interface that allows a direct computation of isochrones 
for any given choice of age, for the chemical compositions and within the age ranges listed 
in tables \M and [2j It also enables the user to compute interpolated evolutionary tracks in 
the mass range 0.5 — 10M & for any chemical composition in the grid. A linear interpola- 
tion on a point-by-point basis between two neighbouring evolutionary tracks stored in the 
database is performed, to determine the track for the specified mass. The mass spacing 
of the tracks included in BaSTI is small enough that a linear interpolation is sufficient to 
guarantee a high accuracy. Once the track is computed, it is immediately transposed in the 
chosen photometric system. 

The use of this interface is particularly simple: first of all one has to get to the server 
using this URL: 

http : //www . oa-teramo . inaf . it/BASTI 

and click on the link Isochrones - Tracks in the Web Tools section; a page like the one 
displayed in Fig. [5] will appear. The user has now to choose among a number of options: 

• the output type, i.e., isochrone or interpolated evolutionary track; 

• the heavy elements mixture, i.e. scaled-solar or a— enhanced; 

• the photometric system of interest; at the moment we allow to choose between the 
standard UBVRIJHKL filters and the ACS-HST ones; 

• core convection treatment, i.e. with or without convective core overshoot; 

• the initial chemical composition (Y,Z). 

Once these input parameters are fixed, the next screen allows the choice of the isochrone 
age (in Myr) or the track initial total mass in solar units, and the job submission. Outputs 
(in the same format of the tracks and isochrones stored in the database) are directly sent to 
the user's browser. 



1 In case of Network and/or technical problem with this official web site, the user can connect to the D. 
Cordier's personal site: http://astro.ensc-rennes.fr where the BaSTI web-tools have been mirrored in 
the BaSTI Web Tools section. 
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4. Luminosity Function Program 

The isochrones stored in the BaSTI database, or obtained with the web-tool described 
in the previous section, can be used as input for the tool described here, that computes 
the luminosity functions (e.g. differential and cumulative star-counts as a function of the 
magnitude in a given wavelength band) of single-age, single-metallicity stellar populations. 

To run this web-tool the user has to provide 4 parameters: the number of isochrones 
whose luminosity functions are needed (maximum 20); the photometric filter (in the UBVRIJKLH 
system); the Initial Mass Function (IMF) exponent; the size of the magnitude bins within 
which the star counts are performed (e.g., 0.15 mag; these bins have to be larger than 0.05 
mag). We adopt for the IMF a form of this type: ^jg- = CM~ a where dN* is the number of 
stars formed with mass between M and M + dM, and C is a normalization constant whose 
value is fixed by imposing that the total number of s tars popu l ating the isochrones is equal 



to 10 . A value of a equal to 2.35 corresponds to the ISalpeterl (119551 ) IMF. 



One has then to choose between the options of computing the luminosity functions only 
until the RGB tip, along the whole isochrone until the Early-AGB, or including also the 
TP-AGB phase. Star counts are then computed analytically by convolving the IMF with 
the distribution of the initial stellar mass values along the isochrones. For example, the 
number of stars N between two consecutive points i — 1 and i along a chosen isochrone is 
given by: N = ^(M}~ a - M}l?). 

As a final step one has to upload the isochrone files, and submit the job. Results 
are directly sent to the user Internet browser as ASCII data, that are divided into several 
sections, each of them corresponding to one individual isochrones. The input parameters are 
recalled: metallicity, IMF exponent, age, etc. The data themselves are displayed in three 
columns: the mean magnitude of the bin, the differential (log iV - where iV denotes the star 
counts per magnitude bin as a function of the magnitude) and cumulative (log N-\ — where 
iV+ denotes the sum of the star counts from the faintest bin to the actual one, as a function 
of their magnitude) luminosity functions. 



5. Synthetic CMD code 

The computation of synthetic CMDs is often required to interpret observations of re- 
solved or unresolved stellar populations. The simplest form of a synthetic CMD is an 
isochrone, that represents the sequence occupied by stars all formed at the same time and 
with the same initial chemical composition. Clearly an isochrone does not directly contain 
information about the number of stars populating the various CMD branches, nor the effect 
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of photometric errors and reddening, and cannot represent composite stellar populations 
made of multiple generations of stars. To include these effects one needs a synthetic CMD 
generator. In the BaSTI website the user can access the SYNTHETIC MAN(ager) code 



that is an evolution of the version briefly described in iPietrinferni et al.l (120041 ) . This code 
computes magnitudes and colors of objects belonging to a synthetic stellar population with 
an arbitrary Star Formation History (SFH). The program employs a grid of isochrones with 
ages between 30 Myr and 14 Gyr (we employ 49 isochrones extended until the end of the 
TP-AGB phase, for each of the 11 metallicities included in our model grid). 

The SYNTHETIC MAN code is based on a Monte-Carlo algorithm, that allows one to 
include in a simple way observation- related effects like reddening/distance/metallicity dis- 
persion and photometric errors, and to identify on a star-by-star basis the various types 
of pulsating stars that may be present in the synthetic population. Synthetic CMDs and 
observation-relate d effects can al so be efficiently determined by means of ana l ytical integra- 
tions, as done in JPolphinl [2002l l It is worth noticing that ISkillman et al.l (120031 ) applied 



three different methods to determine t he SFH of the galaxy IC1613. One of them is the 



analytical technique by iDolphinl (120021 ) . another one is a Monte-Carlo method similar to 
the one we implemented. A comparison of the results with the different techniques shows 
striking agreement. 

The general structure of the SYNTHETIC MAN code is as follows. An SFH file has to 
be specified first. The SFH contains a grid of N ages t{ (the upper limit for N is set to 200) 
increasing with increasing running index i (the present time is assigned t=0) and for each 
age one has to specify a scale factor SFi proportional to the relative number of stars formed 
at that time. Also, the metallicity (parametrized by [Fe/H]) of the stars formed at that 
age has to be given ([Fe/H] J together with a la Gaussian spread around this value. The 
scale factor and [Fe/H] for ages between U <t < are then assumed to be equal to SFi 
and [Fe/H] i; respectively. In other words, the SFH is considered to to be a sequence of step 
functions. The scale factor and [Fe/H] at i = N are not considered and can be set to any 
arbitrary value; the last point in the SFH file is important only because it provides the upper 
age-limit for the simulated stellar population. If one or more single burst stellar populations 
are needed, one or more pairs of identical ages (i.e. £j = t(i+i)) need to be included in the 
SFH file. 

After the SFH is read by the code, the following cycle starts, running along the SFH 
index, from i = to i=N. The age difference between two generic tj and £(j+i) is computed, 
and for each star formed in this time interval (the number of stars formed in a generic time 
interval is determined from the input SFH) a random value of the stellar age ti < t < t(i + 1) 
is drawn with a flat probability distribution, together with a mass M selected randomly 
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according to a user specified IMF (the default mass range goes from 0.1M© to 12OM , 
although a different lower mass limit can be specified by the user) and a random value of 
[Fe/H] assigned according to the value of [Fe/H^ plus the la spread given in the SFH. With 
the three specified values of t, M and [Fe/H] the program interpolates quadratically in age, 
metallicity and then mass among the isochrones in the grid, to determine the star luminosity, 
effective temperature, actual value of the mass (in principle different from the initial value 
because of mass loss) and photometric properties in the U BVRIJHKL system. 

For each generated mass an additional random number determines whether the star is 
member of an unresolved binary system (the percentage of unresolved binaries has to be 
specified as input); i f this is true, the mass of the second component is selected randomly 



following IWoo et al.l (120031 ) and the fluxes of the two unresolved components are properly 
added. 

Once the photometric properties of an object (single or belonging to an unresolved 
binary system) are determined, depth effects of the synthetic stellar population are simu- 
lated probabilistically according to a uniform stellar distribution with a user specified total 
depth (in mag). After depth effects are included, the individual magnitudes are modified by 
adding the effect of extinction (another input parameter; we have employed the exti nction 



ratios from the Asiago Database on Photometric Systems -see iMoro &: Munaril (120001 )- and 
i?y=3.1) and then further perturbed according to a Gaussian distribution, to simulate the 
photometric errors, with a user specified la width. The star [Fe/H] value is also perturbed 
according to a Gaussian error with user specified la width, to mimic the effect of spectro- 
scopic observational errors. As an additional option, the code can search for variable stars, 
according to their location with respect to the boundary of RR Lyrae and Cepheid instabil- 
ity strips; pulsation periods are then determined. Table [3] shows literature references for the 
relevant type of variables. 

Once all stars formed between ages U and t(i+i) are generated, the next interval < 
t < t(i+2)) is considered, and the cicle continues, ending when all stars in the final age bin 
t(N-i) < t < t(N) are generated. 

As a technical comment, if the value of M for a star generated by the Monte-Carlo 
procedure is too large to be still evolving at its age t, or is lower than O.5M , the contri- 
bution to the total mass of stars formed in the population is taken into account, but the 
star photometric properties are not determined. The total mass of the objects formed, the 
total number of the stars whose photometric properties have been determined (in principle 
different from the total number generated for the reasons mentioned above) plus integrated 
magnitudes (in the U BVRIJHKL system) of the synthetic population are also computed. 



-12- 



Like all Monte-Carlo based simulations, our program m akes an exten sive use of random 



numbers. Our random number generator (initially written by lJamesI (119901 )) needs seeds to be 
initialized; these seeds can be either provided by the user or got from http : / /www . random . org 
which derived true random numbers from atmospheric electromagnetic noise. 

The web interface to SYNTHETIC MAN can been found in the BaSTI section called: 
Synthetic Color - Magnitude Diagrams. This software is as simple to use as the others but 
a registration is required before its first use. The user receives a notification by e-mail when 
the computations are completed. In order to get an user identification, one needs to 
contact S. Cassisii or D. CordieJE In the following we summarize the input parameters 
for a run of SYNTHETIC MAN. 



• A scale factor for the SFH. As explained before, the SFH contains the relative weight of 
the star formation episodes in the various age bins. In order to fix the absolute number 
of objects in each stellar generation, a scale factor has to be chosen. The product of 
this scale factor times the relative weight in the input SFH gives the total number of 
stars created in a generic age bin, with mass between 0.1 (or a user specified lower mass 
limit) and 120M Q . If more than 2 10 6 objects are expected to be formed at a given 
time step, the program will stop to highlight the excessive amount of computational 
time needed for the whole simulation. 

• A choice for the mean photometric error. The user can adopt a la constant error that 
it is applied to all stars in the simulation, in all 9 photometric bands. The size of this 
constant error (between 0.0 and 1.0 mag) has to be specified. Alternatively, the user 
can choice to adopt an error varying with the actual star magnitude and/or photometric 
band. The exact values have to be specified by the user in an appropriate input file. 
This photometric error input file contains 18 columns and an arbitrary number N of 
rows (upper limit iV=200). The first 9 columns display N magnitudes M\ (i running 
from 1 to N) in order of increasing values, for the UBVRIJHKL photometric bands 
respectively. The remaining 9 rows are the corresponding la photometric errors in the 
UBVRIJHKL filters, respectively (e.g., ajj, cxb, ay and so on). The choice of the error 
for any individual synthetic star proceeds as follows. If M v is the star magnitude in 
the V band, the program searches for a pair of neighboring tabulated M v magnitudes 
such that M v < My < M^ +1 , and assigns to the star ay the value tabulated for M % y X . 



2 cassisi@oa-teramo. inaf.it 

3 daniel . cordier@ensc-rennes . fr 
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If My is smaller (larger) than M v (My ) the error corresponding to M v (My ) is used. 
The same procedure is repeated for all other photometric bands considered. 

The la mean spectroscopic error (between 0.01 and 1.0 dex). 

A choice for range considered in the stellar mass random extraction. An upper limit 
of 120M Q is always considered. The user can decide to use a different lower mass limit 
(larger than 0.1M Q ) that has then to be specified by the user. 

The total spatial depth of the population (between 0.0 and 10.0 mag). 

The colour excess E(B — V) (between 0.0 and 10.0 mag). 

The fraction of unresolved binaries. 



The m inimum mass ratio for the unresolved binary systems, that enters the IWoo et al. 
( 120031 ) relationship (between 0.0 and 1.0) 

The IMF. The user can adopt a power law of the form M~ x is employed. The exponent 



x has then to be specified. Alternatively, the user is allowed to choice the lKroupa et al. 
jl993h IMF. 

Choice of the isochrone set to use according to the value of an integer index between 
1 and 8. The correspondence between integer value and isochrone set is reported in 
Table HI 

Possibility to search for pulsating variable stars harbored by the synthetic population, 
and compute their periods. An integer value equal to 1 enables the search, a value 
equal to prevents the search. 

Input the desired SFH. A number of pre-specified SFHs for selected galaxies are also 
provided (see Table ED- 

As all simulations based on a Monte-Carlo method, our program needs seeds for random 
numbers generator initialization. The user can choose the values of the seeds, or let the 
program automatically get seeds from the server of entropy http : / /www . random . org. 
We underline that all computations initialized with the same seeds will lead to identical 
outputs. 



After specifying the input parameters, users can start their computations by clicking 
the "Submit" button. The calculation may take some time (from a few minutes to hours, 
mainly depending on the requested number of stars). After the calculation is completed the 
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user receives an e-mail as notification. As the size of the output files can amount to several 
MBytes, they cannot be sent by e-mail and a solution involving a dynamically generated 
web page has been found preferable. In the output web-page four data files are available: 

• Synth_Pop_*_user . in: this ASCII file recalls the parameters used for the calculation. 

• Synth_Pop_*_ . HRD . gz: this gzipped file contains data for the individual stars, in which 
one can find the following quantities: identification number, log(t) (yr), [Fe/H], the 
value of the actual mass, log(L/L ), log(T cff ), M v , (B-V), (U-B), (V-I), (V-R), 
( V — J) , (V — K) , (V — L) , (H — K), the value of the initial mass, the initial mass of the 
unresolved binary companion (if present - for stars without companion this quantity 
is set to 0.0), logarithm of the pulsation period (in days - if the search for variable is 
off or if the star is not pulsating this quantity is set to 99.99), an index denoting the 
type of variable (see Tab. [3] for explanations). This file can be uncompressed with the 
standard GNU software gunzip. As this file can be relatively large, the downloading 
is automatically forced when users click on the link. 

• Synth_Pop_*_ . sf h: contains the Star Formation History data used for the calcula- 
tions. 

• Synth_Pop_*_ . INT.PROPERTIES: file with the integrated UBVRIJHKL absolute mag- 
nitudes produced by the stars evolving in the synthetic population, several selected 
integrated colors, the total number of stars with computed photometric properties, the 
total mass of stars formed (within the mass range specified in the input file) according 
to the specified SFH and a summary of the various types of variable stars found (if the 
search for variables is off the values are all equal to zero). 

The symbol * denotes a string chosen by the WEB tool, to avoid that output files 
corresponding to a specific simulation are overwritten by those related to a different run. 
The string is selected according to the time when the numerical run has been launched. For 
example, if the run is launched at 16.51.34 of May the 3rd of the year 2006, the string will 
be 'May_3_16. 51. 34.2006' 

An example of the output web-page is displayed in Fig. [6j The V — (B — V) Hess's 
diagram taking into account all the synthetic objects drawn for the simulation is also auto- 
matically displayed, stars numbers being encoded with colors. 

The user should be aware that - for obvious reasons - web pages displaying program 
outputs are deleted by the software manager one month after the end of the computation. 
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6. Conclusion 

In this paper we have presented the content of the BaSTI database for stellar population 
synthesis studies. In particular, we have discussed the recent extension of the models until 
the end of the TP-AGB phase, and described three web-tools to compute synthetic CMDs, 
luminosity functions and user-specified isochrones and stellar evolution tracks. 

In the near future the BaSTI database will be developed following two directions: the 
extension of the database itself and the development of an on-line based version of our stellar 
evolution code. 

We plan to extend the mass range spanned by our database, by adding tracks for very low 
mass objects (below 0.5 M ) and very massive stars (above 10 M ), the latter until C- 
burning ignition. 

With a Web-based version of our evolution code, people who are not expert in stellar evolution 
calculations will be able to compute models with their customized set of parameters (mixing 
length value, extension of the overshooting region, metallicity, etc), avoiding inaccurate 
interpolation and/or extrapolation between and/or from existing models. That could also 
allow for more complex tuning like: opacity tables or nuclear rate table switching, change 
of equation of state, etc. We think that this Web-based evolutionary code, as well as the 
BaSTI database could be relev ant tools for t he whole scientific community within the ongoing 



Virtual Observatory" project (]Pasianll2004l ). 
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things. We wish also to acknowledge all people who have already used and/or will be using 
the BaSTI web-tools for their own research, and who will send us their comments/suggestions. 
We warmly thank E. Sandquist for all the useful comments and for pointing out the existence 
of some problems in a preliminary version of our WEB tools. S.C. warmly thanks financial 
support from MIUR and INAF. 



A. BaSTI Web Tools: some technical aspects 

The World Wide Web (WWW), originally conceived for document delivery, has evolved 
into a medium supporting interactive computations. One way for a web server (in our case 
APACHE, see http : / /www . apache . org) to interact with data-generating programs is the well 
known Common Gateway Interface (hereafter CGI). Input parameters and data provided by 
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the user are passed to our three FORTRAN programs through PERL programs. PERL is 
probably the most common CGI scripting language, for which an huge collection of modules 
is available in the Comprehensive Perl Archive Network (see http://www.cpan.org). 
The tracks/isochrones interpolation program and the luminosity function program are sim- 
ply driven with a specific PERL interface, their outputs are directly displayed as an ASCII 
flux which can be saved in the user's hard disk and used by other programs without any 
difficulty 

The CMD generator is probably - at the moment- our more sophisticated WWW-based 
software. When it is invoked, a computation identification string of characters (CISC) is 
assigned. The files are handled and stored with this CISC as part of the file name -so even 
if multiple users from different parts of the world submit simulations at the same time, they 
can all be processed separately. When the computation is completed, an HTML page is 
created on-the-fly, displaying data output files and a pre- visualization made with GNUplot. 
The Perl module Net allows automatic e-mail sending for user notification. 
Figure [7| summarizes the BaSTI Web-Tools general scheme. Some details have been omit- 
ted, like interaction with the plotting sofware GNUplot and random seeds download from 
http://www.random.org. At the moment all calculations are treated after each user's click- 
ing, and consequently several computations can be running at the same time on our machine, 
an occurrence that could seriously affect the computation time. Depending on the success of 
our on-line programs, we could implement as a new feature a more elaborated computations 
managing system that could "queue" the users requests, decreasing individual computation 
times. 
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Table 1: Initial chemical composition of the model grid. 



Z Y [Fe/H] (scaled solar mixture) [Fe/H] (a-enhanced mixture) 



0.0001 0.245 -2.27 -2.62 

0.0003 0.245 -1.79 -2.14 

0.0006 0.246 -1.49 -1.84 

0.0010 0.246 -1.27 -1.62 

0.0020 0.248 -0.96 -1.31 

0.0040 0.251 -0.66 -1.01 

0.0080 0.256 -0.35 -0.70 

0.0100 0.259 -0.25 -0.60 

0.0198 0.273 +0.06 -0.29 

0.0300 0.288 +0.26 -0.09 

0.0400 0.303 +0.40 +0.05 
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Table 2: Summary of the BaSTI model and isochrone database. 



mixture 


scaled-solar 


a-enhanced 


V 


0.2 


0.4 


0.2 


0.4 







0.2 





0.2 





0.2 





0.2 


N° tracks 


20 


20 


40 


20 


20 


20 


40 


20 


M min (M Q ) 


0.5 


1.1 


0.5 


1.1 


0.5 


1.1 


0.5 


1.1 


M max (M e ) 


2.4 


2.4 


10 


10 


2.4 


2.4 


10 


10 


N° isoc. 


63 


44 


54 


44 


63 


44 


54 


44 


Age min (Myr) 


30 


30 


30 


30 


30 


30 


30 


30 


Age max (Gyr) 


19 


9.5 


14.5 


9.5 


19 


9.5 


14.5 


9.5 


Color-T e// 


UBVRIJHKL - ACS HST 
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Table 3: Variable type index and sources used in the population synthesis program. 



Variables Type Index 


Variable Type 


Source for Instability Strip Boundaries 





no variable 




1 


fundamental RR Lyrae 


Marconi et al. (2003): Di Criscienzo et al. (2004) 


2 


first overtone RR Lyrae 


Marconi et al. (2003): Di Criscienzo et al. (2004) 


3 


fundamental anomalous Cepheid 


Marconi et al. (2004) 


4 


first overtone anomalous Cepheid 


Marconi et al. (2004) 


5 


fundamental classical Cepheid 


Bono et al. (2000) 
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Table 4: Index for the iso chrone sets employed by the SYNTHETI C MAN code. 

Index Isochrone set 

1 Scaled solar, no overshooting, 77=0.2 

2 Scaled solar, overshooting, 77=0. 2 

3 Scaled solar, no overshooting, 77=0. 4 

4 Scaled solar, overshooting, 77=0.4 

5 a-enhanced, no overshooting, 77=0. 2 

6 a-enhanced, overshooting, t?=0.2 

7 a-enhanced, no overshooting, 77=0. 4 

8 a-enhanced, overshooting, 77=0.4 
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Table 5: Available Star Formation Histories. 



Stellar Population 



Source 



NGC6822 (global SFH) 
SMC (global SFH) 
LMC (bar field) 
Local disk 
Milky Way bulge 
Sextans A 
LGS3 



Gallart et al. (1996) 
Harris k Zaritskv (2004) 
Holtzman et al. (1999) 
Rocha-Pinto et al. (2000) 
Molla et al. ( 20001 
Dolphin et al. (2003) 
Miller et al. (2001) 
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-8 - 



-7 - 



-6 - 



-5 - 



-4 - 




Z = 0.01 
Z = 0.008 
Z = 0.004 
Z = 0.0006 



77 = 0.2 



J I L 



77 = 0.4 



J I L 



8.0 



9.0 

log(t) yr 



10.0 



Fig. 1. — Run of the K-band SBF magnitudes of Magellanic Cloud superclusters as a func- 
tion of their age. Theoretical SBF magnitudes for the appropriate metallicity and age range 
(denoted with different line styles) of the individual superclusters (superclusters of different 
metallicities are denoted with different symbols) are also displayed. These theoretical values 
are obtained from the scaled-solar isochrones including overshooting and with 77=0.2, ex- 
tended along the TP-AGB phase. At Z=0.0006 we display also values obtained with r]=0A 
(dashed-dotted line - see text for details). 
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Fig. 2. — As in Fig. [T], but for the J-band SBF magnitudes (see text for details). 



-26- 




8.0 9.0 10.0 

log(t) yr 



Fig. 3. — As in Fig. (U but for the integrated (H — K) and (J — K) colors of Magellanic 
Cloud superclusters (see text for details). 
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Fig. 4. — Cumulative K — (J — K) CMDs of AGB stars in a sample of LMC clusters. 
Theoretical isochrones for the appropriate age and metallicity range are also displayed (see 
text for details). 
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Fichier Edition Affichage Aller a Marque-pages Qutils Aide 

- - 1 r-j) http://ast i-o.ensc-rennes.fr/basti/iso_tracks/index, ht rr 



J ©ok [__r~ 



Musique . __; Librairies -Editions CBanque Olnformatique jMoteurs £3 Linux lz3 Astro _] Dan [j CPI 2 Qj MailGaelle [_ LJ1300 Q CUPS 



WEB TOOL I 



U... 1 1. Ttvtranlnr 


Cordie 


r from ENSCR 


France! ii 



from Ljvorpq 
PMetnnrV-n'i] .- 



is a web interface written in F5KL language by Et. [' 
collaboration with Dr. S. Cassisi from Teramo 
'istrinferni from Tn~?i> ■:''j.>^;^|-.-rii:...rv (Italy) and Dr. M. Salaris 
'Jiiiv^r-vii' v (UK). This program uses evolutionary tracks from 
alan:. M. and 'Ostein F. i2Q04. Ap.I -vol. 612. 158) . This Web 



Tool is not a simple applet and as it is server side based, it is robust and fast. 




'BVRUHKL [Scaled solar or Alpha enhanced) 
CS [Scaled solar only] 



■ Standard [without overshooting) 
Non standard (with overshooting) 



= 0.0001 Y= 0.245 
Z= 0.0003 Y= 0.245 
Z= 0.0006 Y= 0.246 
/ 0.001 V 0.246 
Z= 0.002 Y 0.246 
Z= 0.004 Y= 0.251 




Fig. 5. — Web interface for the tool dedicated to isochrones and evolutionary tracks. 
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The BASTI Population Synthesis Program Web Interface - Mozilla Firefox 



Fichier Edition Affichage Ajler a Marque-pages Outils Aide 



E3 



| ^ http :ffastro ,ensc-rennes ,fr/basti/synth_pnp/result5/Synth_Pnp May_15_l 1 ,49,25_: " | © OK | 



tJMusique- DLibrairies -Editions i_JBanque ■ J Informatique OMoteurs LiJLi'nux IZJAstro U Dan U CPI 2 U MailGaelle u LJ13D0 fj CUPS 



aSTI 



outputs "Synth Pop_May 151 1.49,25 2006 " foi 



daniel.eordier-AT-ens c-rennes .ft 



• The parameters file: Synth Fop May 15 11 4R25 2005 us~r.m 

• The Gzipped synthetic HRD file: Svnth Pop May 15 11.49.25 2005 .HRD.qz ( size: 312KB) 

• The stellar formation history file used for this calculation: Svnth Pop May 15 11.49.25 2006 ,sfh 

• The integrated colors file: Synth Pop May 15 ll.4fi.25 20C"? .INT COLORS 



BaSTi - Population Synthesis CHD 



interface created by D. Cordier 




[Hess's diagram of the output, the scale on the right gives the number of stars per pixel.] 



Termine 



Fig. 6. — Example of the output web-page for the SYNTHETIC MAN web-tool. 
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BaSTI user 




http : //www. oa-teramo . inaf . it/BASTI 



Fig. 7. — Data fluxes in the BaSTI web tools. Random number generator seeds provided by 
http : / /www . random . org has been omitted. 



